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Using a metric-based formalism to treat cosmological perturbations, we discuss the connection 
between anisotropic stress on the brane and brane bending. First we discuss gauge-transformations, 
and draw our attention to gauges, in which the brane-positions remain unperturbed. We provide a 
unique gauge, where perturbations both on the brane and in the bulk can be treated with generality. 
For vanishing anisotropic stresses on the brane, this gauge reduces to the generalized longitudinal 
gauge. We further comment on the gravitational interaction between the branes and the bulk. 
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One outstanding problem in brane world cosmology is to develop a better understanding on the evolution of 

Q ■ perturbations. To gain insights is of prime importance. In fact, in order to make predictions regarding the primordial 
spectrum of perturbations or the anisotropies in the cosmic microwave background, the creation and evolution of 
£C) i perturbations have to be understood. Only then we can hope to test the theory with future cosmological experiments. 
On the other hand, the question of stability of the global brane world space-time has to be addressed, because some 
(or most) solutions of Einstein equations could be in principle unstable. To attack these issues, the use of a full 
five-dimensional description is necessary. 

Some work has already addressed fluctuations in brane world theories. Namely, creation of perturbations on 
the brane was considered in Q-Q, where specific assumptions were made from the beginning. Other groups have 
developed a formalism to treat perturbations in brane worlds for rather general situations, see [|J (hereafter Paper 
£Nj ■ I) and 0] In Paper I, a metric-based formalism, which represents a straightforward extension of the usual 

four-dimensional approach Jl3| ], was developed for a brane world theory with two branes. Also, some of the evolution 
equations on the brane were found. However, two assumptions where made, namely: i) it was required that the 
position of the branes was not affected by a first order perturbation, and ii) a generalized longitudinal gauge (GLG) 
was used in the bulk. A result of these assumptions was that the anisotropic stresses on the branes had necessarily 
^ 1 to vanish. Although, this particular scenario would be adequate to study scalar fields or ordinary matter without 
anisotropic stresses, it does not represent the most general case. The aim of this paper is precisely to fill this gap by 
considering the most general case of having anisotropic stresses on the branes and by allowing the brane positions to 
be perturbed. 

We will consider only the case of single extra dimension, which is assumed to be compactified on a circle S\ with 
a ^-symmetry. The (unperturbed) branes will be located at the fix points of the ^-symmetry, which, in our case 
are taken to be y\ = and y2 = R- 

The paper is organized as follows: in the next section we will discuss gauge transformations. After briefly reviewing 
the construction of gauge invariant variables, we will show that the GLG is not compatible with unperturbed brane 
positions in the most general case. We will then discuss a novel gauge, which reduces to the GLG for vanishing 
anisotropic stresses on the branes. We will also discuss the Randall-Sundrum gauge in the context of cosmological 
perturbations. In Section III, we will derive the perturbed Einstein equations for the gauge introduced in Section II 
and we will derive the junction conditions. We will also briefly discuss the gravitational coupling between the branes 
and the bulk in Section IV. Finally, our conclusions can be found in Section V. 
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II. GAUGE TRANSFORMATIONS AND GAUGE-INVARIANT VARIABLES 



In this section we begin with a brief review of the gauge-invariant formalism for metric perturbations in brane- 
world models (see Paper I). Then we will discuss gauge-transformations and the inclusion of the anisotropic stresses. 
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Throughout this paper, we will consider the case for five-dimensional brane worlds only. The additional spatial 
coordinate is denoted by y 5 = y. For the purpose of this section, all we need to specify is that the branes are 
stretched across the usual four-dimensional space-time and that they are located at specific points along the additional 
dimension. We will be more precise below. 

The most general higher-dimensional metric consistent with the maximally symmetric three-dimensional spatial 
manifold is given by (here the indices a and b are either or 5 and y = t) 

ds 2 = a 2 { lab dy a dy b - <>.,,/.,•' r/.r' } , (!) 

where the scale factor a and the metric "f ab are functions of the coordinates y a only and is the metric of the 
three-dimensional maximally symmetric space, given by 



1 + - x l x m S h 



(2) 



with k = 0, 1, —1 for flat, closed or hyperbolic 3-geometries, respectively. Given this structure of the background, we 
are able to classify metric perturbations by their three-dimensional tensor properties as in the four-dimensional case 
fl3f . The perturbed metric can be generically written in the form 

ds 2 = a 2 { lac (S c b + 2<f> c b ) dy a dy b - [(1 - 2i\>) fly + 2E {ij + 2F m + fcy] dx l dx 3 - 2W al dy a dx 1 } , (3) 

where Fi and hij have a vanishing divergence and hij is traceless. Also, the three-vectors W a i can be split as follows, 

W ai = B a]l + S m , (4) 

where S^u = 0. 

An infinitesimal coordinate transformation can be written as 

x a -> x~ a = x a + (5) 

with (infinitesimal) parameters £ Q . We split these parameters as = (£ a ,£ 4 ), where £ % ' = rf + ^ l being ^ a scalar 
function and if a divergenceless three-dimensional vector, i.e. tfu = 0- We use the convention that indices of type a 
(i) are lowered, raised and contracted using the metric jab i^-ij)- Furthermore, we take the vertical bar to denote the 
covariant derivative with respect to "f ao or Oy depending on the index type. As shown in the transformation of 
the scalar sector is found to be 

Hob = -S,(a\b) ~ H c £ c j ab) (6) 

H = H a ta, (7) 

6E = (8) 
5B a = £ a -Z {a . (9) 

where we have introduced the generalized Hubble parameters 

H c = ^. (10) 
a 

The vector perturbations in the metric ((^) change according to 

SFt = -Vi, (11) 
SS ai = -ru\ a ■ (12) 

The tensor perturbation hij is, of course, invariant under the first order gauge transformation. From these transfor- 
mation rules one can easily construct gauge invariant variables for scalars and vectors, which we repeat here (see 
for further details): 

Scalar variables 

$ab - <f>ab + H C {B C - E\ c ) lab + {B {a - E\ (a )\ b) (13) 
* = 4> - H c (B c - E\ c ) . (14) 
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Vector variables 



In this paper we shall consider scalar perturbations only, for which the perturbed five-dimensional metric (||) reduces 

to 

ds 2 =a 2 {b 2 [(1 + 2(f>)dt 2 - 2Wdtdy - (1 - 2T)dy 2 ] 

- - 2tp) + 2E\ ij ] dx l dx J - 2B a \ l dtdx l - 2B 5{i dydx i } (16) 

where we have defined 

= ^8, r = -4, w = 2 ( f> 5 = -2cf> 5 . (17) 

It is convenient to use a conformal gauge for the metric 7 a b, i.e., 

( 7o6 ) =6 2 diag(l,-l). (18) 

Note that, as it is further explained in Paper I, this gauge for 7 a fc can always be chosen. Here b = b(t,y) is a new, 
independent scale factor. 

So far the discussion was independent of the existence of branes. We now introduce brane sources, which are, in 
general, located at specific points in the fifth dimension. As in Paper I, we are going to consider the setup motivated 
by heterotic M-theory. The extra dimension is compactified on the orbifold Sx/Z^. That is, we compactify the fifth 
dimension on a circle restricting the corresponding coordinate y to the range y G [-R, R] with the endpoints identified. 
The action of the Zi symmetry on the circle is taken to be y — > —y. Associated to this symmetry there are two fixed 
points at y = y\ = and y = j/2 = R, where we assume that two three-branes, stretching across 3 + 1-dimensional 
space, are located. 

We have to truncate the five-dimensional metric in order to make it consistent with the orbifold symmetry. For the 
unperturbed metric this implies the following conditions at the fix points: 



9nu{-V) =9ixu(y), (19) 

g^{-y) = -9v&{y), (20) 
g&si-y) = 555(2/)- (21) 

When the background geometry is perturbed, however, the branes are in general no longer located at y n =constant, 



but at a perturbed positions which depend on the intrinsic brane coordinates x % and t. Thus, the symmetry conditions 
(|l9|)-(^T|) cannot directly be applied to the perturbed background. However, by a suitable coordinate change, one 
can always gauge away these perturbations on the brane positions (see the discussion below). Therefore, using these 
particular coordinates, (fl9|)- (pil ) also hold for the perturbed metric. 

Note that the conditions (|19|)-(]2l|) impose some restrictions on the gauge transformation (j^) that we are allowed to 
perform, if the new coordinate system is assumed to have the same symmetries. Indeed, in this case we have to make 
sure that gauge transformations do not lead out of the class of metrics defined this way. In particular, this means 
that the parameter £ Q for an infinitesimal coordinate transformation has to be restricted by 

e{-v) = e{y) (22) 
e(-y) = -e(y), (23) 

which directly follows from and the symmetry conditions ([l9|)-(^l|). ^From these rules we can deduce the 

Z2 properties of the various scalar quantities in metric ([l6]): while the background scale factors a, b as well as the 
perturbations 0, ip, T, E and Bq are Z2 even, that is, for example, a{—y) = a(y), the perturbations W and B5 are 
Z2 odd, that is for example W(— y) = —W(y). Similarly, for the scalar components in the transformation parameter 
£ Q , we find that £0 and £ are even while £5 is odd. 

To be precise, as we did in Paper I, we are going to work in the boundary picture, where instead of working with 
the whole orbifold, only a half of the circle is considered. We will require that all components of metric (nq) are 
continuous across the full orbicircle except for the odd components W and B$ which are allowed to jump at the fixed 
points (but continuous otherwise). To simplify the notation, we define the value of an odd field on the brane as the 
one that is approached from within the interval y G [0, R\. This is precisely the boundary value of the field as viewed 
in the boundary picture and it represents one half of the jump of this field at the fixed point. 
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A. Relocating the brane positions: a gauge for non— bending branes 



In certain cases, the formalism of brane world perturbations is much simpler if the position of the branes remain 
fixed at constant y. Indeed, if the branes do not bend into the fifth dimension, then the Zi symmetry along the 
orbifold will be straightforwardly imposed by the conditions (|l9|)-(pTl). Furthermore, as explained in Paper I, the low 
energy-effective action may easier be obtained if the brane positions are unperturbed. 

In this subsection we will show that it is always possible to gauge away the perturbation of the brane position. 
We then motivate a novel gauge, in which brane world perturbations can be treated very generically. We will start 
considering a gauge for a background with two branes located at fixed y = y n positions. Then, the most general first 
order perturbation in this scenario will be described by the metric (^|), together with two first order perturbations 
on the brane positions. That is, the branes will bend into the fifth dimension according to, 

y ( h n) =y n + X n (t,x i ). (24) 

Hence, the most general set of first order perturbations is <f>, if), T, W, E, B a and X n . In that case, the ^-symmetry 



applies at the new brane positions y b and, thus, the simple conditions ( 19 )— (|2l|) no longer apply. However, this set 



of perturbations is too large. In fact, observe for instance that we can construct two new gauge invariant quantities 
out of Ei 5 , B 5 and t/£ ( n= L2) by the following simple combination, 

B 5 -E l5 ~y { b n \ (25) 

which directly follow from (||)-© an d the transformation rule for yj,: Syb = £5 after a first order gauge transformation 
(|J). Thus, if we choose B§ — Ei 5 = everywhere in the bulk, i.e. if we choose the GLG, both brane positions are in 
general perturbed. However, as it will become clear in the following section, only when the anisotropic stresses on the 

branes vanish, can we use the GLG and yet keep the branes at y b = const.. Hence, in the GLG the position y^ b n \x^) 
of the branes contain only information about the anisotropic stresses. 

Observe that, without loosing essential degrees of freedom, it is always possible to perform a gauge transformation 
that resets the perturbed branes to their initial unperturbed positions y — y n . The procedure is the following: 

The general metric (|l^) contains five degrees of freedom associated to the arbitrary choice of the five quantities 
that describe a five-dimensional gauge transformation. However, in order to relocate the brane positions, we need only 
to partially fix the arbitrary function £ 5 . In fact, since under a general gauge transformation £ a the fifth coordinate 
y changes according to y = y + £ 5 , we can take, 

e(x a )=e^ a )+i 5 (* Q ), (26) 

being £| (x a ) any odd function across the unperturbed brane positions, but smooth otherwise, with the following 
boundary conditions at the unperturbed brane positions 

£(t,x\y n ) = T K(t,x t ), (27) 

where the upper (lower) sign applies when the position y = y n is approached to the right (left). On the other hand, 
^ 5 (x a ) is any smooth function with the following boundary conditions at the unperturbed brane positions, 

i 5 (t,x i ,y n ) = 0. (28) 

Note that the boundary condition ( p7|) allows to gauge away the bending of the brane by fixing again the position of 
the brane at y — y n (see Fig 1.), while the term £ 5 with boundary ( p8| ) contains a residual gauge freedom which still 
remains after fixing the position of the branes. Obviously, any further gauge transformation should satisfy (|28|). The 
reason is because the ^-symmetry applies again to the unperturbed brane positions y n . Thus, £ 5 is an odd function 
across y = y n and to avoid the appearance of distributional terms on the metric coefficients, £ 5 must indeed be zero 
at y = y n . 

While we relocate the position of the branes, we may choose two smooth functions £ = E and £0 = E\ — Bo by 
means of which E = Bo = 0. Then, we have the following generic set of perturbations (f>, ip, T, W, E — 0, Bo = 
and by (|) 

B 5 = B 5 + £5 - £| 5 = B 5 - E ]5 + £a 5 + | 5 - (29) 

However it is not in general possible to choose a quantity £5, that vanishes on both branes, by means of which B5 = 0. 
Indeed, to do so we must choose £5 = Ei 5 — B5 — ^5. However, the odd quantities Ei 5 — B 5 — £\ 5 do not in general 
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vanish on the branes. Therefore, the combinations ±(E^' — + \ n /b 2 ), which represent half of the jump of the odd 
functions — _B 5 — £\ 5 across the branes, are already gauge invariant. Indeed, we have not further gauge freedom 
to change them. Obviously, this is a simple consequence of (Eq) being gauge invariant. On the other hand, as we will 
explicitly see in the next section, E}^' — -Bg — £a 5 are directly related to the anisotropic stresses on the branes. 




FIG. 1. In the most general situation, the brane positions along the fifth dimension y are perturbed (a). However without 
affecting the essential gauge freedom, it is always possible to gauge away the perturbation on the brane positions. Once the 
branes have been relocated, we can use the remaining gauge freedom to set E — Bo — and simplify B$ as much as possible. 
Observe that, by relocating the brane positions at y = constant, we can straightforwardly impose the Z2 symmetry of the 
orbifold and then the junction conditions for the metric coefficients can be read off directly from Einstein equations. 

Although it is not possible in general to switch to a gauge where the odd perturbation B$ is zero, we can actually 
eliminate most of it as follows. Since it is always possible to choose B$ in the form (|2^), we can then split 

B 5 =B 5 -E l5 + L 5 + £ 5 , (30) 



5 



where L5 is the broken line that connects the values of B5 on both branes (from within the interval y £ [yi, 2/2]) and 
shares the same jumps of the function B§ across the branes. On the other hand B$ — E\ 5 is a smooth function which 
vanishes on the branes and which determines the separation between the functions B§ and L§ (see Fig. 1). The 
broken line L5 can be easily constructed as follows 




(i)+ 



B 



(2)- 



B 



(1)+' 



2/2 - yi 



(y - vi) 



{9(y - yi ) - 9(y - y 2 )} 



where the plus (minus) sign in holds for the value of B% when the brane is approached to the right (left). Now, 

explicitly using that B§ 



(n)± 



±{B^ ] - E\"> - A„/6 2 ), such a broken line can be easily simplified as follows 



1. — (flf-^-^l+i; 



6 2 



R (n) p (n) 
-°5 ^15 



B 



(n) 



^|5 ^2" 



y-yi 
yi - yi 



(31) 



Observe that all the physical information that the perturbation B$ carries is actually encoded in L5. This is true 
because, B§ — is & smooth function that vanishes on the branes and thus it can be gauged away by the simple 
gauge transformation 



& = E {5 - B 



(32) 



in (£5(j). Then, after this last gauge transformation, the value of the perturbation S5, which for simplicity we will 
denote by B, would be reduced to its simplest value, namely 



B 



(33) 



We note here, that this particular gauge reduces to the usual GLG when the jumps of the metric perturbation B§ 
vanish on the branes. We will see below that this is the case when the anisotropic stresses on the branes vanish. 

With all this information in mind, we should be careful with the construction of the complete set of gauge invariant 
variables discussed in the previous section. Of course, in order to avoid distributional terms in the construction of the 
gauge invariant variables, then (|l3|) should be constructed using B$ — E\ 5 instead of B5 — E\$. Indeed, that would 
prevent the appearance of delta-terms coming from the y-derivative of the odd quantity B$ — E\ 5 , which does not in 
general vanish on the brane positions. 

Observe that, since B§ — _B| 5 is taken to be zero in this gauge, the metric perturbations <f>, ip, T, W actually 
correspond to the gauge invariant variables. On the other hand, B is determined by the broken line (|3lj), whose slope 



depends on the variables Bg — E^ 1 ' , which are defined on the branes only and are gauge-invariant under coordinate 
transformations which respect the symmetry conditions (|l9|) — (pl|). 

In summary, by means of a gauge transformation compatible with leaving the branes at their unperturbed positions, 
we have been able to eliminate two degrees of freedom (E and Bq) in the metric ( |l6| ) instead of the three degrees 
of freedom (E, Bq and B§) that can be cancelled with an ordinary gauge transformation. The remaining degrees of 
freedom that cannot be cancelled are the jumps of B$ — E\ 5 + £>5 across the branes, which, as we will see in the next 
section, are directly related to the anisotropic stresses on the branes. 

For simplicity in our notation we will denote the set of gauge invariant variables in the gauge discussed above by 0, 
-0, r, W and B. It is worth noticing that, since it is always possible to keep the two branes at fixed y — y n orbifold 
positions, such a gauge provide great advantages in solving the full five-dimensional dynamics. As it will become clear 
in the following section, the dynamics of the variable B is described by the evolution of the anisotropic stresses on the 
branes only. Thus, when B = 0, this particular gauge will reduce to the generalized longitudinal gauge introduced in 
Paper I. 
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FIG. 2. In general the metric perturbation B5 cannot be gauged away by means of a gauge transformation which is compatible 
with the Z2 symmetry along the orbifold. The reason is that B5 is an odd quantity and thus it may jump on the brane positions 
(a). Since the Z2 symmetry conditions ([l9|)-(|2"l|) forbid the jumps of B$ from being gauged away, these jumps are already gauge 
invariant. However, any behavior of the perturbation B$ in the bulk besides the purely linear is totally irrelevant. In fact, it 
is always possible to perform a gauge transformation compatible with the Z2 symmetry which reduces B5 to the straight line 
L5 in (a). This is simply possible by choosing a gauge transformation with £5 = — (B5 — L5). Observe that such a £5 vanishes 
on both branes (b) and thus is compatible with the Z2 symmetry. 



B. Other possible gauges 

In the previous subsection we have used the gauge freedom to cancel the perturbations E and Bq and to simplify 
the perturbation B5 as much as possible and yet keep the branes at their unperturbed positions. We have also said 
that this particular gauge choice corresponds to the generalization of the four-dimensional longitudinal gauge. On 
the other hand, another choice for the gauge fixation that is commonly found in the literature is setting the five- 
dimensional perturbations T, W and B 5 to zero, and keep the perturbations E and B instead, i.e. the the scalar 
part of the Randall-Sundrum gauge. However, this choice of gauge is generally incompatible with keeping the branes 
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at their unperturbed positions. In fact, let us suppose that we start from the most general set of perturbations, i.e <fi, 
ip, r, W , E, Bq and -B5, and we want to set T to zero. In order to do that, we will use the transformation equation 
(pfi, which can be written as, 



b 2 ST = £' 




(34) 



Observe that, regardless of the particular value of the gauge fixing function if we want to set the perturbation T to 
zero we need to solve a differential equation for £5 with the boundary conditions £5™^ = on both branes. However, 
this is not in general possible. As a particular example, consider the Randall-Sundrum case where 6=1 and a = a(y), 
then equation (|3~i| ) would simply be, 



with general solution, 



&-~&-r = o, 

a 



C(t,x l ) + / T(t,x\y) 



dy 



(35) 



(36) 



where C(t, x l ) is an arbitrary integration constant. Observe, thus, that it is always possible to fix C(t, x l ) to set £5 = 
in one of the branes but not in the two of them at the same time. Thus, suppose that we fix C(t,x l ) in order to set 
^ 2 - ) = on the right brane located at a constant y = y%. Then on the left brane £W 7^ 0, which means that this 
brane cannot be located at a constant y — y\ position but, according to (|^), its new position should necessarily be 



y = Vi-^5{t,x\ yi ) (37) 

where £s(i, x i , y±) is given by evaluating (|3^) at y = y\. Then, any gauge transformation (|) from this new position 
has its fifth component zero in agreement with the ^-symmetry. A similar situation would happen if we tried to set 
the perturbation W = 0. 

In summary, although it is possible to use the Randall-Sundrum gauge for brane world scenarios with a single brane 
located at a fixed y position, it is not in general possible in this particular gauge to fix the positions of more than one 
brane. Such an effect was already pointed out in |ll| . 



III. THE FORM OF THE PERTURBED EINSTEIN EQUATION AND JUNCTION CONDITIONS 

In this section, we calculate Einstein equations in the gauge considered in section 2. The metric in this gauge is 
given by 

ds 2 = a 2 {b 2 [(1 + 2(f>)dt 2 ~ 2Wdtdy - (1 - 2T)dy 2 ] 

- - 2tp)} dx l dx 3 - 2B\ i dydx i } . (38) 

Note that the only difference between such a gauge and the generalized longitudinal gauge considered in Paper I, is 
the quantity B — B§. As mentioned in the last section, and explicitly seen below, B = B§ contains only information 
about the anisotropic stresses on the branes. 

The two branes arc located in the orbifold at y = constant and the Einstein equations can be written as 

1 2 

G af ) = R a p - -g a pR = T a p + ^2 T 5i )( % ~ Vn), (39) 

1 n=l 

where we have set the five-dimensional Newton constant to one, for simplicity. The delta-functions in this equation are 
covariant with respect to the fifth dimension, that is, they include a factor of 1/^—555. Furthermore, T a p is the bulk 
stress-energy tensor induced by fields that propagate in the full five-dimensional space time. The brane stress-energy 

tensors T^l , on the other hand, originate from fields that are confined to the branes at the orbifold fix points. We 
note here, that we are not considering Gaussian normal coordinates. 
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Next, we specify the stress-energy tensors. As described in Paper I, they can be found using standard methods. 
For the bulk, the most general form of the unperturbed stress-energy tensor is 




and the background stress-energy tensor on the branes has the form 

,(«) o 

T(n)a 
P = 



,(n)jt 





pv'O'j 




For the perturbed stress-energy tensor we find 

dp 

ST a = ' 



rB\* + (p + p)v^ 



-(p + p)b 2 v b 
-8p5 l 3 - 



Sr + 2r((j) + r) -{p + q)W {p-q)b- 2 B\ 



13 

b~\ 



\.i 




(40) 



(41) 



(42) 



(43) 



where v and u are two potentials for "velocity" fields and a, satisfying ct' 1 ^ = 0, determines the anisotropic stress. 
The perturbed brane stress-energy tensors are 

/ 6pM -( P W +pW)b- 2 v\f -SrW 

§T (n)a /3= ( p („) +p (n)} v (n)\i _ Sp (n) § i. + a {n)\i _ u {n)\i 

^ Sr^-p^W p^b^By-b^Jf 

where, unlikely in our previous work pL we have explicitly included an anisotropic stress term the brane 

stress-energy tensors and a non-vanishing component 8T^ 1 ^. The reason is that, as we will see below, once we 
include a perturbation B in our formalism we get new delta terms which should be precisely matched with a brane 
anisotropic stress and a brane 5T^ 5 component. We would like to present the equations of motion based on the 
metric ( p8| ) and on the above stress-energy tensors that follow from the Einstein equation (^). The background 
equations following from ( |39| ) have already been given in ref. 0] and thus we should refer to this work for details. 
However, since the introduction of a B perturbation will also introduce new terms on the perturbation equations at 
linear order, which are not present in our previous work B, we explicitly give them here: 
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(44) 
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(45) 
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(46) 







-(p + p)v -J2(P (n) +P (n) )v {n) ~8{y - y n ) , 
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(47) 
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3 d d b d Id 2 
2adt + bdt~ 2dt? 



(q-p)B + u + J2 (-p (n) B + U W) S(y -y n )\ , 



(48) 



(a6) 2 <5G l 



a' 9 „d 9 6 9 6' 9 d 2 



a b b 2 a dy a dt bdt b dy dt 2 



d b d b 1 d d 2 



b'b V 

2 2 6 3 3 

b 2 b a a dt a dy b dy b dt dtdy 



2 2^ 6^ 3°'^ 3^^ ^ ' ^ ^ 

b 2 b a a dy a dt b dt b dy dy 2 



W 



+ 
+ 



a 2 b 2 



6^A_6-— 2— -2— 2kb 2 
a dy adt dy 2 dt 2 



V + 2& 2 (^ - + r) 



. a 



d 

a ay 

:> - + f 

a ay 



B^-b 2 ^-^ + rf y 



z 



(n)l* 



<%- Vn) 



(49) 



Here we have defined the delta- function S which incorporates a factor (—555) = 1/ab. 
We write down also the traceless part of eq. (691): 



n 1 ! i -i»-* + r)V' i = -«V |j -^(3^ + |- 



\j 3 



'j-" Ife 



(50) 



Observe that, for the particular case where the anisotropic stress in the bulk vanishes, i.e. cr' 1 ^ = 0, we may identify 

+ -* +r — p^ + sj)* (51) 



As a next step we analyze the background and the perturbed Einstein equations: 

Background equations 

Matching first the delta terms in the background equations we can easily get the junction conditions for the scale 
factors a(t, y) and b(t, y) (see Paper I for details): 



a 6 2 V 



(52) 



Here and in the following the upper sign holds for the brane at y = y\ and the lower sign for the brane at y = 1/2 • 
On the other hand, from the 05 and 55 component of the background equations, we find 
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p(») = -3- + p (n) ) T 2a6r , (53) 



2 7,2 



a ab , , o a b 

T + kb 2 = 

a ab 



1 p {n) (p {n) + 3pW ) ! y 



12 



(54) 



which respectively represent the background energy conservation equation on the brane and a background dynamical 
equation for the scale factors on the brane. 

Perturbed equations 



Observe first that the components SG° , 6G l i , SG 1 j (i ^ j), SG° i , SG° 5 and SG 4 contain explicit delta-function terms 
and they should be matched by terms containing first y derivatives of W and B and second y derivatives of all other 
quantities. This leads to 

^' = -W-l (4 n) - O ' " ± \<* (sp {n) - rp (n) ) , (55) 



a + b + Ot r ± 3 



ia6 (jp^ - ry™)) ± iafe (V n > - r P ("') , (56) 

±y ° {n) \ = - - + ^ K» - < } ) " fe , (57) 

W = -i (flW - T J + P (n) ) (58) 

5 r (n) 

^=-7^- (59) 



p(") 



I* 



(n) 

Tt 



In these equations the background junction conditions ( p2| ) have been used. From eqns. ( p8|) and (|59| ) we find 

(nl 

5r („) = T « p( „) ^(»J + ^ „(„) _ P_ ^(n) _ 

Recall that the quantities <5r( ra ) and are uniquely fixed by the other components and they are, generally, non- 
vanishing, as already discussed in Paper I. In fact, they are directly related to the anisotropic stresses on the branes. 
For instance, we can easily transform (61) into an equation involving the anisotropic stresses on the brane as follows: 

= (ofcrWy T JpW (pW (62) 

where we have introduced for convenience the traceless quantities SIL^ j and y(*0 j, which are given by 

5^\ = 5r^%- 5 hr^\, (63) 

V,(nr. =v (n)lj._| v (n)lf fc< (64) 

Finally, we would like to discuss the influence of the new terms discussed in this section (i.e. the quantities B = B$ 
and cr'™-') for the low-energy effective theory discussed in Paper I. Because B in the bulk simply connects the values 
of the anisotropic stresses on the branes (see the line-construction in section 2), the y— averaged B along the bulk is 
zero. However, at the position on the branes, the jumps of the function have to be taken into account, so that the 
effective four-dimensional anisotropic stress 04 has to be replaced by 



n=l 



where < ... > describes the averaging over the fifth direction. The other terms are not affected (R is the coordinate 
distance between the branes and % is a modulus describing the size of the fifth dimension) . 
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IV. GRAVITATIONAL INTERACTION BETWEEN THE BRANES AND THE BULK 



In this section, we would like to discuss the gravitational interaction between the branes and the bulk. Since our 
gauge is well adapted to analyze the whole five-dimensional brane world, we should be able to get some insights on 
the non-local nature of this particular geometry. 



Observe, that the variables Sr^ are the components 8T^°5 of the stress-energy tensor on the branes. Thus, 



ii.-, 




is connected to an energy flux, i.e. it determines a flux of energy along the fifth dimension. 



Observe from fl63| ) that the anisotropic part of Sr^ is directly related to quantities like (p^ n ' +p^)V^ n ' > or <jw, which 
are naturally related to the third order time-derivative of the reduced quadrupole moment of the matter fields on the 
brane^. Observe, that the emission of gravitational radiation by a cosmological source is directly linked to the third 
order time-derivative of its reduced quadrupole moment (see |25|| for details). Therefore, we may physically interpret 
the anisotropic part of the quantities Sr^ as a non-local measure of the energy exchange between the branes and 
the bulk through the exchange of five-dimensional gravitational radiation. 

The generalization of the energy-momentum conservation equation ( p3| ) for brane perturbations can be found by 
restricting the 05 and 5i-components of the perturbed Einstein equations to the branes. This respectively gives the 
continuity equation and Euler's equation. These two equations together with (|6l]) give the full energy-momentum 
conservation for the brane matter field perturbations under possible interactions with the bulk. In fact, making 
explicit use of the junctions equations these three equations may be conveniently written as follows 



IfipW + 3-(*pW + 5p ^) - 3 V n) + P (n) )V; - -(P 
at a a 



(n) 



P 



± 2ab(p + q)6r {n) + 2ab(Sr + 20r + Tr) , (66) 



0_ 

Ot 



2 2- 



=F 



ab 



(67) 



aV(p - q)Bf + b 2 kB^ + ^B^" ]ik + \aH\ 



0_ 

dt 



= T ab (pW+p(. n A v W 



(68) 



These equations provide a coupled system of differential equations for the matter perturbations on the brane. In 
particular observe that: 

i) Equation (|6q ) describes the growth of pressure and density perturbations sourced by r, Sr, Sr^ and a Poisson- 
like term v , fc from the velocity perturbations. 

ii) Equation ( |67| ) describes the growth of velocity perturbations, which are sourced by 5p( n \ (f>, the anisotropic 
stresses on the branes and the momentum flux from the bulk matter. 

iii) Equation (|68| ) describes the growth of anisotropic stresses on the branes sourced by SA™" 1 and velocity pertur- 
bations on the branes. This equation states that anisotropic stresses can be smoothed out by changing into 
velocity perturbations and gravitational radiation. For the particular case that Sr^ — 0, then equation ( |6~l| ) 
states that decays on the anisotropic stresses can induce velocity fluctuations on the matter fields only (reference 
[ p2[ provides details on this type of effects for ordinary cosmology). 

Let us collect the physical meaning of the quantities r, Sr and Sr^ n ': 

i) r determines a non-perturbed energy flow in the bulk arising from the 05-component of the background stress- 
energy tensor (EG). For the particular case that the bulk is empty, i.e. p + q = 0, this term is zero. 



1 The reduced quadrupole moment of the matter fields on the brane is defined as, P = f p \x % x^ — ^S lj 
is related to the non-spherical part of the kinetic energy of the source which is related to the anisotropic stress. 



x k x k ] 



d J x. Thus, r 
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ii) 8r determines a first order perturbation on the energy flow in the bulk arising from the 05-component of the 
perturbed stress-energy tensor (0). Also, for the particular case that the bulk is empty this term is zero. 
Recall, however, that being Sr = does not mean that there is not a flow of energy through the bulk. In 
fact, it is possible to have pure gravitational radiation propagating in the bulk and yet have Sr = 0. This is 
true because Sr comes from a purely local stress-energy tensor for first order bulk perturbations. However, 
the energy associated to gravitational waves is non-local. Indeed, first order gravitational waves are vacuum 
solutions of Einstein equations and their energy can be only properly computed by non-local space-time averages 
of quadratic wave amplitudes (see for instance )0. 

iii) 8r^ determines a first order perturbation of the "projected" energy flow on the branes, which may have two 
contributions: a local one (due to the matter distribution) and a non-local one (due to gravitational waves). 
While the local projection of Sr^ on the brane is zero when the bulk is empty, in general SA n ^ ^ due to the 
non-local contributions. It is worth noticing that for domain walls we have 5r n = 0, because the equation of 
state is p + p = and the anisotropic stress vanishes. This suggest, that, at least for the scalar part discussed 
here, domain walls do not radiate energy (at first order), see also p3[. 



Finally, the restriction of the 55 component of the perturbed Einstein equations, eq. (45), to the brane leads to the 
following evolution equation 



6 2 (2^-0)'V- 3^-3-0 + 3 ( 7-3-1 ij + 6kb 2 (i' + (^) + a 2 b 2 p {n]2 [i(l + 3u; (n) ) 
1 a \ b a I [b 

Sq + 2ql If 3 {r A s(n) S^_ s{n) ar 



. (69) 



V. DISCUSSION 



Using the formalism presented in Paper I, we have discussed how anisotropic stresses on the branes can be easily 
included. We have also argued, that if we use the generalized longitudinal gauge in the bulk, the branes no longer 
remain at fixed y = const, orbifold positions. On the other hand, the brane positions are related to the anisotropic 
stresses only. Also, we have introduced a novel gauge, in which the positions of the branes may be fixed at y — const, 
while the information on the anisotropic stresses is encoded into the metric perturbation B. Finally, we have discussed 
the gravitational interaction between the branes, which, in our formalism, is described by the 05-component of the 
brane energy momentum tensor. 

An outstanding problem is to find solutions of the bulk Einstein equations (Q)-(^). It is worth noticing that some 
work has already addressed this issue from the point of view of a brane-observer only. However, in the presence of 
perturbations, this is not enough. In fact, one needs to solve the full five-dimensional Einstein equations in order 
to fully understand the non-local interaction between the two branes. Some work along this direction was already 
done in Paper I for the linear M-theory background. A derivation of the low-effective perturbation equations in 
Randall-Sundrum brane world [^6| (with two branes) is still lacking. However, even without specifying the details of 
the theory, one is able to derive rather general statements about the evolution of perturbations and the consequences 
of the dynamics of the perturbations, (see Q, Jj7|, || Q, |15| and for a discussion on scalar perturbations, |l^] 
for a discussion on vector perturbations and pl | for tensor perturbations). For example, in fli] ] it was shown that 
due to the energy-flow, onto or away from the branes, will modif y t he evolution of super-horizon amplitudes on the 
brane. In fact, in contrast to the usual results in four dimensions [p~3f| , adiabatic perturbations will not be constant in 
general. As discussed in Jl4[ , this effect is a consequence of the statement that super- horizon amplitudes are constant 
if the energy-conservation holds |2jJ) . However, energy-conservation can easily be violated in brane worlds. 

Large scale anisotropies for a brane in an Anti-de Sitter bulk were discussed in |l9"| . Also, using the gauge presented 
in this paper, the full set of equations on the brane for this scenario was derived in pO[). Hence, there is the hope, 



2 The propagation of linear gravitational waves is determined by the first order perturbation of the vacuum Einstein equations. 
However, the backreaction of the gravitational waves over the background geometry is determined by Einstein equations with 
a non-zero stress-energy tensor associated to the waves. Such a tensor can be obtained by non-local averages of quadratic 
combinations of the first order gravitational wave amplitudes. 
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that the formalism presented here can be used to fully understand the interaction between the bulk gravitational field 
and the brane matter. 

Note added: While this paper was prepared for publication, ref. ]Tl} | appeared, which addressed some of the 
questions in section 2. 
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